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Abstract. We find the minimum scale factor, for which the nonnega- 
tive Bottcher-Wenzel biquadratic form becomes a sum of squares (sos). 
To this we give the primal and dual solutions for the underlying semi- 
definite program. Moreover, for special matrix classes (tridiagonal, back- 
ward tridiagonal and cyclic Hankel matrices) we show that the above form 
is sos. Finally, we conjecture sos representability for Toeplitz matrices. 



1 Introduction 

The Bottcher-Wenzel inequality ([1], [2], [3], [7], [6]) states that for real square 
matrices P, Q the biquadratic form 

BW = 2 (||P|| 2 ||Q|| 2 - trace 2 (P T Q)) - ||PQ - QP|| 2 (1) 

is nonnegative, with the Frobenius norm used. Replacing the factor 2 by 2+y n , 
it is natural to ask for the minimum y n such that 

(2 + Yn ) (llP|l 2 HQII 2 " trace 2 (P T Q)) - ||PQ - QP|| 2 

is a sum of (polynomial) squares. We answer this question in Theorem 1 by 
showing that the minimum value is (n — 2)/2. 
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For simplicity, we use one-subscript notation for the entries of P and Q, 
described by means of the "small" index matrix 

IND = ((i-l)n + j)ft =1 . 

Then P and Q can be generated by vectors p and q of dimension m = n 2 as 

P(i,j) =p(IND(i,j)), 1 <i< j <n. 

Introducing an index matrix will be especially useful in Sections 3 to 5, where 
tridiagonal, backward tridiagonal, cyclic Hankel and general Toeplitz matrices 
will be investigated. For these special cases we prove (for Toeplitz matrices: 
conjecture) that the corresponding BW form is a sum of squares (sos). 

It is quite odd that although (real) Hankel matrices are symmetric, thus 
normal, hence nonnegativity easily follows [1], this does not imply that a sos 
form also exists (except if n = 3, Example 9). On the other hand, Toeplitz 
matrices are usually not normal, yet the corresponding BW form is sos, at 
least according to our well-grounded Conjecture 15 at the end of the paper. 

2 The case of general matrices 

Let P, Q be arbitrary n x n real matrices with entries 

P = (P(i-1)n+j)iJ> Q = ( £ l(i-1)n+j)y> 

as indicated above. (Notice that we use this indexing technique for simplicity.) 
It turns out [5] that the above forms depend only on the variables 

zy =Vi<\) ~ qtPj> 1 < i < j < tv, 
i.e. on the skew symmetric matrix 

Z = pq T -qp T 

of order n 2 , a benefit of including the term trace 2 (P T QJ. Indeed, we have 

||Z|| 2 =2(||P|| 2 ||Q|| 2 -trace 2 (P T Q)), 
and all entries of the commutator [P, Q] obviously are linear forms of the Zn-S. 
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Let us formulate the primal and dual semi-definite programming problems 
(see e.g. in [8]) for the eigenvalue optimization: 

min {tr(CX) : X > 0, tr(AtX) =0, 1 < i < M, tr(X) = 1} (Primal) 

M 



max{-y M +i: S = C - Y y t A t - u M+ i I > 0} 



(Dual) 



t=l 



where C, S, X, A t and the identity I = In are all real symmetric Nth order 
matrices, C and (A t )^ are given, the primal matrix X, the dual (slack) matrix 
S and vector y are the solutions of the program, tr(AB) = trace(AB) denotes 
the scalar product of the symmetric matrices A and B, and > stands for the 
semi-definite ordering: A > B iff A — B is positive semi-definite. 

The quantities zy will play the role of 'candidate monomials' (better to say, 
differences, and hereafter called candidates) with ordering 

Z=(Z1,2> 21,3, Z2,3> Zl,4j ■■■> Z\,n, ■ ■ ■ > Zn-l,n) T - 

The indices can be read from the "big" index matrix 



/0 1 
• 



POS 



2 

3 




4 7 

5 8 

6 9 
10 



\ 



of order n 2 to be 



k = i + 



(j-i)0-2) 



1 < t < j < n 2 



Note that zy = for all i, and that the entries below the diagonal are omitted 
due to zy = — Zj^, enabling us to reduce the number of unknowns. Also note 
that IND is related to P and Q, while POS is connected with Z. 

As an example, we give the biquadratic form BW as a quadratic form of the 



quantities (zm) for n = 3. Observe that ||Z| 



1<i<j<9 z 'i,j- 



Example 1. For n = 3 the objective takes the form 



BW = ||Z|| 2 -(Z2,4+Z3 )7 ) 2 -(z 1) 2+Z2,5+Z3 )8 ) 2 -(Zi ) 3 + Z2,6+Z3 )9 ) 

— (— Z1,4 — Z4 ; 5 + Zsj) — (— Z2,4 + Zg^) 2 — (— Z3,4 + Z5,6 + Z6,?) 2 

- (-Z1,7 - Z 4 ,8 - Z 7 ,9) 2 - (-22,7 ~ ^5,8 ~ Zs,?) 2 ~ (-^3,7 + Z6,s) 2 - 
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(Note that the form BW can be thought of as a function of the matrices 
(P, Q), the vectors (p, q), the matrix Z, or of the vector z.) We give now all 
the (quadratic) relations holding for the variables (zy)i<i<j< n 2 as 

Zy Zk,i + zt,i Zj.ic - Zi,k Zj,! = 0, 1 < I < j < k < I < n 2 . (2) 

These easily checked relations define M = (^ ) symmetric constraint matrices 
At, each having exactly 6 nonzero (off-diagonal) entries. For instance, equation 

^2,3 ^4,5 + Z2 > 5 Z 3,4 ~~ z 2,4 ^3,5 = 

defines an A t with nonzero entries in positions (3, 10), (8,6), (5,9) and their 
transposes, see the matrix POS. Now we can state our main theorem. 

Theorem 1 The minimum value o/y n , for which (1) is a sum of squares is 

n-2 

Yn = — ^— • 

Proof. We give the optimal primal and dual solutions and describe the main 
characteristics of the optimal dual matrix. Since the objectives coincide, the 
strong duality theorem yields the desired result. 

By fixing the order of the variables (zy) above, matrix C is uniquely de- 
termined. To get the (slack) matrix S = C — J~ UtAt, we use the following 
strategy. Note that we not only give the set (At) of active constraints (as e.g. 
when taking the half Newton-polytope), but also give their coefficients (yt)- 

Strategy A. Assume the commutator [P, Q] contains an entry (zy +Zv \+. . . ) 
with i, j, k, I distinct and i < j, k < I. Then the quadratic form z T Cz associated 
with C necessarily contains a term 2zuZ\^\. We 'halve' this term, and leave 
one Zy Zv i unchanged as is, while apply for the other term the basic quadratic 
relation (2). By using the correct sign, this defines a constraint A t and the 
corresponding dual variable tjt for some t. Finally, let yra+i = 

Now we give the obtained primal and dual solutions. In view of the quite 
combinatorial character of the problem, we do not detail each block, instead 
we give some explanations and important cross-references (control sums) and 
for matrix S we provide Table 1. with all essential informations. 

The Primal Problem 

Before defining the optimal primal matrix X, we note that its rank is Q). 
For indices (i, j) : 1 < i < j < u we define the vectors Vy of dimension ) 
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to have 4(n — 1 ) nonzero coordinates (four 2's and 4(n — 2) ±1 's) using rowi, 
rowj, column^ and columnj of the index matrix IND, cf. Example 2 below. 
Next we form the matrix of these vectors 

V=[v 1)2 , v 1;3i v 2)3 , v n _i, n ], 

and define matrix Xo = V V T = Y_ v ij v J j with the following properties: 

Xo is a symmetric matrix of order N = ) and rank Q) . The ■ 's are 
orthogonal with norm square ||vy || 2 = 4 • 4 + 4(n — 2) • 1 = 4(n + 2). The trace 
of Xo is 

tr(Xo) = X tr(vy v?j) = IK if = 4 ( n + 2 ) (*) ' 

i<j i<j ^ ' 

thus by defining 



n 



X=(4(n + 2)ri] ] X 

we get a trace 1 matrix. The eigenvalues of Xo are 4(n + 2), those of X are 
(2) 1 (hence Q) ] X is a projection). The Vy's are also eigenvectors of C: 

Cv y = (2-n)/2 Vy. 

Furthermore we have 



|Vij 



|2 



tr(CXo) = ^trfCvyvJ,) = ^vJjCvy = L 

i<j t<j i<j 

and finally, the primal objective equals 

f(fY1 tr(CX ) 2-n 

tr(CX) = = — r— . 3 

tr(X J 2 

The DuaJ Problem 

The matrix S = C — 2ZytAt resulting from Strategy A is positive semi- 
definite and decomposes into some blocks given in Table 1. (Observe that all 
eigenvalues and diagonal entries of 2S are integer - a reason for the factor 2.) 

Here we list the important facts and control sums concerning the blocks of S, 
and in Example (3) we give further hints for understanding the construction. 

ROW-control: an element in the last column of row i is the scalar product 
of row 1 and row i. For instance, the number of zero eigenvalues — the defect 
of S— equals to Q)2 + 1 (n - 1 ) = n 2 - 1 . 
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1 


No of blocks 


(2) 




3(4) 


(2) 


Total 


2 


Block sizes 


6n-8 


(?) 


4 


1 




3 


Eig=0 


2 


n-1 






u 2 -1 


4 


Eig=4 


1 








n(n-1)/2 


5 


Eig=n 


2n-4 


(V) 


1 




(n 2 -1)(n-2)(n + 4)/8 


6 


Eig=n + 2 


3n-5 




2 


2 


(Ti/-1)(n-2)(n + 4)/4 


7 


Eig=n + 4 


n-3 




1 




n(n- 1)(n 2 +n-2)/8 


8 


Eig=2n + 2 


1 








n(n-2)/2 


9 


Diag 


n(+2) 


n-2 


n + 2 


n + 2 


(u J -n)(u z +2n-4)/2 



Table 1: Decomposition of the matrix "2S" 



EIG-control: the last column (the number of eigenvalues, rows 3 to 8) sums 
up to ), the order of the matrix 2S. 

DIAG-control: The sum of the elementwise products of row 1, 2 and 9, 

n \ (2n*n+4(n-2)(n + 2))+lQ (n-2) +3 (^\{n + 2) + Q 1 (n+2) 

equals to (n + 1 ) (n 2 + 2n — 4) , the trace of 2S . (In the blocks of order 6n — 8 
there are 2n diagonal elements "n", and 4(n— 2) diagonal elements "(n+2)".) 
TRACE-control: the trace of the coefficient matrix C equals 

tr(C) =2( 11 " 2 ) -ufn-1)-fu 2 -n)u = nfu-1) 2 fn + 1 



The first subtrahend comes from the diagonal of the commutator [P, Q] , the 
second from their off-diagonal elements. Due to diag(S) = diag(C) +y n I, the 
connection between the traces of matrices C and S is 



tr(2S)=2(tr(C) + ^(^f) 



The number of all constraints is (4 ) , while that of active constraints equals 

n ( n 2 ') + n{n ~ " ( ( n 2 2 ) + 2{n ~ 2) ) = (2) lnI ~ 41 s 3(n + 21 (3) ■ 



Here the first term is associated with the main diagonal of R = [P, Q] (by 
virtue of Z\a = there are only n — 1 terms in R(i, i)), while the rest comes 
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from the off-diagonal of the commutator R (where there always are two terms 
for which the basic relations do not apply, see the Example 2). 

Thus S is positive semi-definite with defect n 2 — 1 , and its eigenvalues range 
in the interval [0, tl + 1]. To sum up, the primal objective (3) coincides with 
the dual objective yjyi+i, the negative of y n , which proves the theorem. □ 

There holds no strict complementarity, for rank(X) = Q) < n 2 — 1 = def(S). 

Example 2 To define the primal matrix X take the four scalar products 

(rowi, colj), (colt, rowj), (rowt, rowj), (colt, colj) 

in the index matrix IND where each of the four products determine n coordi- 
nates in vy as follows. If n = 3 and i= 1, j = 2, then rowi = [1,2,3], C0I2 = 
[2,5, 8] T which yields by (1,2) ~ 1, (2,5) ~ 8, (3,8) ~ 24 the coordinates 
1, 8, 24, see also matrix POS. Similarly we calculate the other three triples, 
giving together 

1, 8, 24; 4, 10, 21 (!); 4, 8, 13; 1, 10, 28. 

The repeated elements (1, 4, 8, 10) denote positions with value 2. The excla- 
mation sign refers to an entry —1 (since (7, 6) must be inverted to (6, 7) ~ 21 ). 
To sum up, we get 

vi, 2 = (2, 0, 0, 2, 0, 0, 0, 2, 0, 2, 0, 0, 1 , 0, 0, 0, 0, 0, 
0, 0, -1 , 0, 0, 1 , 0, 0, 0, 1 , 0, 0, 0, 0, 0, 0, 0, 0) T . 

Example 3 Hints for obtaining the dual matrix. We give some details for the 
(2) most important blocks of order 6n— 8. There is a one-to-one correspondence 
between these blocks and ordered pairs (i, j), 1 < i < j < n. To collect the 
indices for the block containing zy, we have to consider the 4(n — 1) terms in 

(rowt, colj), (colt,rowj), (row^rowj), (colt, colj) 

(the same as for vy above!) and further 2(n — 2) terms in the products 

IND(iJ) * diag(/ i, j), IND(j,i) * diag(/ I, j). 

Here diag(^ i, j) stands for the n — 2 entries of the diagonal of IND, differing 
from i, j. As in Example 2, choosing n = 3, i = 1 , j = 2, vector diag{^ 1,2) 
reduces to the (3,3) entry 9, thus we get (using index matrices IND and POS) 

((1, 2), (3,3)) ~ (2,9) -30 and ((2, 1), (3,3)) ~ (4,9) ~ 32. 
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Hence the diagonal block containing row 1 (related to Z\z) also contains rows 
30 and 32. The whole index set at issue is [1, 4, 8, 10, 13, 21, 24, 28, 30, 32], 
and the corresponding block is the 10x10 (irreducible) matrix 



( 


3 





-2 





-1 





-1 








\ 







3 





-2 





1 





-1 










-2 





3 











-1 


-1 













-2 





3 


-1 


1 
















-1 








-1 


5 








-1 


-1 


1 







1 





1 





3 


-1 













-1 





-1 








-1 


3 
















-1 


-1 





-1 








5 


1 


-1 
















-1 








1 


5 





V 














1 








-1 





5 / 



with eigenvalues (0, 0, 3, 3, 4, 5, 5, 5, 5, 8). 



3 Tridiagonal (and backward tridiagonal) matrices 

In a former paper [4] we have shown that for nth order matrices P, Q with only 
nonzero entries in row 1 and column n the BW form is sos, however in case of 
(additional) main diagonal elements this is no more true. Therefore one would 
guess that 3n + 0(1) nonzero elements cannot be allowed, however the result 
below shows that the answer depends on the position of these elements. 
We shall use an index matrix given e.g. for n = 3 as IND =(345 



.0 6 7. 

Lemma 4 For tridiagonal P, Q the BW form is sos, especially we have 

bw = 



- y~ (Z3i-4,3i-3 - Z3i-l,3i) 2 - X. z 3i-2,3i-1 _ X. Z 3i,3i+3 

_ X U3i-2,3i-l + ^3i-1 ,3i+l ) 2 — X (z 3 t_2,3i + Z 3i . i3 i + i ) 2 

= y~ (Z3i-4,31-1 + Z3i-3,3i) 2 + (z 3 i_ 4;3i — Z 3i _ 3)3 i._i ) 2 

+ X (z 3 i_2,3i-l — z 3i-l ,3i+l ) 2 + ( z 3i-2,3i — z 3i,3i+l ) 2 

+ X z 3i,3i+2 + ^ X Z 3i _ 2)3 i+1 +2 X Z tj~^_ Z 3i-l,3i+3- 

i+5<j 

Remark 5 The first equality gives the biquadratic form at issue, the second 
one is the claim: the sum of squares representation. (The negative terms in 
the last row are evidently canceled.) 
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n 


A = 


A = 1 


A = 2 


A = 3 


2-bl. 


act. 


rk(X) 


2 


3 





3 





2 





3 


3 


7 


1 


12 


1 


6 


1 


7 


4 


11 


2 


30 


2 


10 


2 


11 


5 


15 


3 


57 


3 


14 


3 


15 


6 


19 


4 


93 


4 


18 


4 


19 


7 


23 


5 


138 


5 


22 


5 


23 


8 


27 


6 


192 


6 


26 


6 


27 



Table 2: "Tridiagonal matrices 1 



Although SDP is not needed here, for the identity of the Lemma can be proved 
directly, we yet give some facts. The eigenvalues of the dual matrix S for the 
actual semidefinite programming problem are integers (0, 1 , 2, 3) in this case, 
too. This is so because matrix S decomposes into at most second order blocks 
of the form (] j ) and ( ^ ^ ). 

Table 2 illustrates the main features of the underlying semidefinite program. 
First the number of the eigenvalues of S are given, then the number of 2 x 2 
blocks in S (the number of scalar blocks is not shown), the number of the active 
(y t 7^ 0) constraints, and finally, the rank of X. (The n — 2 active constraints 
correspond to the positions (i — 1 , i), (i, i — 1 ), (i, i + 1 ) and (i+ 1,i) in IND.) 

It is easy to get a formula for these quantities, e.g. the number of eigenvalues 
At = 2 can be determined by subtracting the number of all other eigenvalues 
from the order (3n - 2)(3n - 3)/2 of S. The result is 3 + 9( n ~ 1 ). 

Note that strict complementarity does hold: the number of zero eigenvalues 
of S coincides with rank(X), the number of nonzero eigenvalues of X. 

Backward tridiagonal matrices 

They have many similar properties, except that the case n odd is worse: 
while for n even all the eigenvalues of S are integers (lying in [0, 4] ) , for n odd 
this does not hold, therefore we write '— ' instead. Also, in this case there are 
(apart from the scalar and 2x2 blocks) 4x4 blocks, too. All this information 
is contained in Table 3 from where one can see that for n even we again have 
strict complementarity, as in the tridiagonal case. 
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n 


A = 


A = 1 


A = 2 


A = 3 


A = 4 


act 


2-bl. 


4-bl. 


rk(X) 
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3 





3 








2 








3 


3 


8 










3 


3 


7 


5 


4 


13 


4 


25 





3 


6 


3 


8 


13 


5 


20 










11 


6 


17 


18 


6 


25 


6 


81 


2 


6 


14 


6 


18 


25 


7 


32 










30 


9 


27 


30 


8 


37 


8 


173 


4 


9 


22 


9 


28 


37 



Table 3: "Backward tridiagonal matrices 1 



Example 6 We calculate the number act of active constraints: 

{5n — 12, n even 
5n — 8, n odd 

The number of terms in a typical row of the product of backward tridiagonal 
matrices usually equals (0, . . . , 0, 1 , 2, 3, 2, 1 , 0, . . . , 0). However, in case of the 
commutator PQ — QP there are some minor changes: for odd order 1, for even 
order 2 main diagonal entries contain only two terms (instead of 3), due to the 
identity zy + Zy = 0. On the other hand, if n is even, there are two opposite 
entries (with indices (k, k + 1) and (k + l,k), where k = tl/2) which do not 
generate any constraint, for the corresponding indices are not distinct. 

Now we easily calculate the number asked, which is e.g. for n = 6 equal to 
5u — 12 = 18. To this consider the matrix 

(1 2 1 0\ 

2 3 2 1 

12 2*10 

1*221 

1 2 3 2 

\0 1 2 2/ 

with the number of terms in a given position of [P, Q] , and take Q) for any 
entry e > 1 . They sum up to 2* (6(1) + (I)) = 18. The general case is similar. 
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4 Cyclic Hankel matrices 

When investigating Hankel matrices, we find that - except for the case n = 3, 
see below - they do not generate sos BW forms. However, cyclic ones behave 
well. We make use of the small index matrix (given for n = 3): IND =^23?^. 

Theorem 7 For cyclic Hankel matrices P, Q the BW form is a sum of squares. 

Proof. Using the above-defined index matrix with (1,2, ...,n) as first row 
and (n, 1 , . . . , n — 1 ) T as last column, we obtain 

I|P|| 2 =^IIP|| 2 , IIQH 2 = n||q|| 2 , trace(P T Q) =np T q, 
consequently 

2 f||P|| 2 ||Q|| 2 - trace 2 (P T Qj) = 2n 2 (||p|| 2 ||q|| 2 - (p T q) 2 ). 



The commutator [P, Q] is a skew symmetric cyclic Toeplitz matrix having 

k = k r " n " 1 ' 



different entries t^ = t[ n ' with row one as 

(0, ti , . . . , t k , -tic, . . • , -ti ) (n odd ) 
(0, ti , . . . , tic, 0, -tic, . . . , -ti ) (n even). 

Thus the subtrahend is ||R|| 2 = 2u^ t 2 , and the whole BW form equals 
Observe now that all terms in 

n— i i 
ti = t t n Zj,i+j - Zj,n-i+j 

J=1 j=1 

are distinct (i = 1 , . . . , k), hence the Cauchy-Schwarz inequality in conjunction 
with nk = n^2yij < (£) imply 

k k n— i i 

Z t ^L n (Z z f,i + i + Z z f,n-i + j) <n!4, 

i=1 i=l j=l j=l i<j 

which proves the theorem. The last inequality turns into equality for n odd. 
□ 
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Remark 8 The case n = 4 is especially interesting. Then the commutator is 

I t -t \ 
-t t 
-t t 

V t o -t o / 

with t = ti = t| 4 ' = zi 2 + 22,3 + Z3,4 -21,4, therefore the formula 

4 (Zl,2 + zf (3 + Z2, 3 + zf,4 + ^2,4 + ^3,4) = 
+ (Zl ,2 + Z2,3 — Zi ,4 + Z3,4 ) + Ul ,2 ~ Z2,3 + Z] ,4 + Z3,4 ) 
+ (Z1,2 + Z1,3 - Z2,4 - Z 3 ,4) 2 + (z 1>2 - Z 1)3 + Z 2 , 4 ~ Z 3;4 ) 2 
+ (Z1,3 + Z2,3 + Zi,4 + Z2,4) 2 + (Z1,3 — Z2,3 — Zi,4 + Z2,4) 2 , 

(a consequence of Eulers identity) yields the sos-representation needed. 

Example 9 The case of (general) third order Hankel matrices. The index 
matrix IND is now y 2 3 4 J , the order of C, S and the constraint matrices {A t } 
is (2) = 10, the number of the A t -s is (A = 5. By help of vector 

Z= (Z1,2, Z1,3, Z 2) 3, Zi,4, Z2,4, Z 3;4 , Zi ;5 , Z2,5) Z 3)5 , Z 4)5 ) T 

and matrix C the objective can be written as BW = z T Cz = tr(Czz T ), which 
becomes - by means of a standard SDP relaxation - trace(CX). Our MATLAB 
program yields^ = (0,0, 1,0,0,0), i.e. only one constraint will be active, giving 



S = C-y 3 A 3 = 
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-1 
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{-1} 











-1 
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-1 
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{-1} 
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-1 
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{1} 





-1 








-1 











1 ) 



(In the original C the six entries in braces are zero.) The last zero iny indicates 
the sos representability. To obtain the concrete sos form, we calculated the 
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eig en- decomposition of the three blocks 



( 


1 


-1 


-1 


1 \ 




-1 


4 


-2 


-1 




-1 


-2 


4 


-1 


V 


1 


-1 


-1 


1 / 



B, 



with integer eigenvalues 

E, : (0 4 6 ) , 



,B 2 



/ 


2 


-1 





-1 \ 




-1 


3 


1 


-1 







1 


1 





V 


-1 


-1 





2 J 



,B 3 



E 2 : (0 1 3 4), 



and (unnormalized, integer, columnwise) eigenvectors 



V! : 



/2 
1 
1 

\0 



1 



1 



1 

-1 
o / 



v 2 



/ 1 



V 1 



i 


2 
1 



1 




-1 



1 \ 

-3 
-1 

1 / 



V, 



2 
-1 



'1 3 



We sum up the result: with zn = p^q j — q^pj, 1 < i < j < 5 the following 
identity holds for the B W form generated by two third order Hankel matrices: 



2z\ 2 + 3zf j3 + 6z| 3 + 2z 2 j4 + 4z| 4 + 6z 



■f z 2 c + 2z 2 



"(Zl,3 +Z2,4 +Z3,5) 



+ 3z| 5 +2z| )5 



3,4 T^l^T ^2,5 
(Z1 ,2 + 22,3 + Z3,4) - (Z2,3 + Z3,4 + 2:4,5) - 
Ul,2 — 22,3 — Z3,4 + Z4,5) 2 + 3(Z2,3 — Z34) 

1 3 1 

+ gUl,3 + 2z 1)5 + z 3i5 ) 2 + -(z 1)3 - z 3;5 ) 2 + -(z 1)3 - 3z 2)4 - z 1)5 + z 3i5 

1 , ,,3 



+ -(Z1 )4 +Z 2> 5) +2^1 ,4 



Z2,5j 



5 Toeplitz matrices 

In this section P and Q will be arbitrary real Toeplitz matrices. Observe that 
the main diagonal entries don't play any role (to prove this use temporarily 

the more standard notation P = (pj-0 and Q = ( q j t ) , then the entry in 

R = [P, Q] equals Y. z k-i,j-k> while Zoj_i + Zj-i,o = 0)- Hence we can reduce the 

number of variables to get e.g. for n = 3 the index matrix IND = ( 3 1 V 

\4 3 0/ 

Another speciality is that now there occur repeated terms as well. To handle 
these, introduce the multiplicity vector p of dimension m by 

Pi = {the number of occurrences of p^ in P, 1 < i < m}. 

Then the following easily proved representation holds. 
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Lemma 10 



m— 1 m 



2(||P|| 2 ||Q|| 2 -trace 2 (P T Q)) =l£_ f_ ^Wi,j> 



i=i j=i+l 



and (since the commutator is skew persymmetric) , ||R|| 2 = 2^ i=1 H 




In view of the lemma, we define the symmetric matrix C by help of equation 



strategy of choosing the dual variables. 

Strategy B. Since the entries of [P, Q] are linear forms in [z\\ ) , their squares 
figuring in ||R|| 2 involve some mixed products of the form ±2zijZici. Whenever 
finding such a term with distinct {i, j,k, I}, we increase the actual value of y. 

It turns out that Strategy B works for n, 3 < n < 7, however, for n > 8 the 
dual matrix S = C — ^_yt^t will have (more and more) negative eigenvalues. 

Lemma 11 For orders n not exceeding 7, the matrix S is p.s.d, i.e. for these 
values the BW form is sos. Some further properties of S of arbitrary order n 
are: the minimum off-diagonal entry of S is — L^y^-J; the defect of S, i.e. the 
multiplicity of zero as eigenvalue is n— 1; the maximal diagonal entry and also 
the maximal eigenvalue is n(n — 2). Moreover, S is a direct sum of two types 
of submatrices of the following order: 

- type (a): 2, 4, 6, . . . ,2(n — 2); (denote by B the largest block here) 

- type (b): 1 , 1 , 2, 2, . . . , n - 2, n - 2, n - 1 . 

The orders of these matrices sum up to (n — 1)(2n — 3), the order ofS. 



with D diagonal, H Hermitian (i.e. symmetric), both of order n — 2. The 
diagonal elements of D are (i(i + 1)) in reverse order: ((n — 2) (n — 1 ),..., 6, 2). 
Matrix H is also of a special structure: the elements on the border are — 1, 
those on the 'neighboring' border are —2, and so on. This matrix is p.s.d. for 
n < 7, but has at least one negative eigenvalue for n > 8. 

Remark 12 The further submatrices of type (a) also are critical, e.g. the next 
one (of order 2(u — 3)) has a similar form with diagonal elements (i(i + 2)) in 
D, while H is the same (of the appropriate size). Therefore there is a second 
negative eigenvalue for n, 14 < n < 20, and so on. In general, the symmetric 
matrices H are of the same form, and the diagonal entries o/D are (i(i + k))^. 




The largest block B of type (a) is crucial. It has a decomposition B 




Sum of squares representation 
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Example 13 Matrices of order 5. In this case P and Q have m = 2(n— 1) = 8 
nonzero elements, the candidate vector z has dimension (™) = 28, the total 
number of constraints is (™) = 70, and the number of active constraints is 14. 

It always suffices to examine the first row and the first column of R, for all 
other entries are contained in these, e.g. R(1 , 1 ) = zi ; 5 + Z2,6 + Z3,7 + Z/^s, and 
R(2, 2) = Z2,6 + 2.3,7. The number of the active constraints for n = 5, coming 
from row 1 and column 1 is indeed 6 + 2 (3 + 1) = 14, as stated above. This 
can be proved by induction, by noting that 




As regards they coordinates, since the product 2zz^5,7 occurs two times (as 
the above formulae show), we write —2 in the suitable positions (overwriting 
the — 1-sJ to get S(13, 17) = S(3,21) = —2, and so on. 



Now we give another example illustrating the role of the crucial block B . 

Example 14 The case u = 8. The matrices D and H are now: 

1 -1 -1 -1\ 

2 -2 -2 -1 

3 -3 -2 -1 
3 -3 -2 -1 " 
2 -2 -2 -1 
1 -1 -1 -\) 

The characteristic polynomial of the block B = factorizes into piP2> 

where w[x) = x 6 -100x 5 +536x 4 -53472x 3 +327472x 2 - 575680* -145 152, and 

pi has a negative zero: -0.2228. (All other zeroes o/pi and P2 are positive.) 

Finally we mention that although the above strategy works only up to n = 7, 
the standard semidefinite program yields results indicating that BW can be 
sos for some larger orders, too, hence we guess that BW is sos in general. The 
difficulty is that the corresponding dual variables (yt) of the program are not 
recognizable real numbers. Nevertheless we formulate the following. 

Conjecture 15 The Bottcher-Wenzel form (1) generated by two real Toeplitz 
matrices is sos, i.e. a sum of squares of polynomials, now: quadratic forms. 
Give - if possible - a rational certification, i.e. rational parameters (yt) such 
that S = C — Y_ ytAt is positive semidefinite. 



D 
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6 








2 ) 



, H 



-1 
-2 
-2 
-2 
-2 
-1 
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